Charge Resistance in a Majorana RC circuit 
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We investigate the dynamical charge response in a "Majorana Coulomb box" realized by two 
Majorana bound states hosted at the ends of a mesoscopic topological superconductor. One side 
of the wire is coupled to a normal lead and low frequency gate voltage is applied to the system. 
There is no dc-current; the system can be considered as an RC quantum circuit. We calculate 
the effective capacitance and charge relaxation resistance. The latter is in agreement with the 
Korringa-Shiba formula where, however, the charge relaxation resistance is equal to h/2e 2 . This 
value corresponds to the strong Coulomb blockade limit described by a resonant model formulated 
by Fu [PRL 104, 056402 (2010)]. We also performed direct calculations using the latter model and 
defined its parameters by direct comparison with our perturbation theory results. 

PACS numbers: 73.43.-f, 74.45.+C, 73.23.-b, 71.10. Pm. 



Introduction. — The search for Majorana fermions 
(MFs) and Majorana bound states (MBSs) in condensed 
matter systems is currently under way, driven both by 
their unique physical properties and by their poten- 
tial applications in topological quantum computation 
P, [H • There are several promising systems predicted to 
host MFs, including: spin-triplet p-wave superconduc- 
tors, topological insulator/superconductor hetcrostruc- 
tures and semiconductor/superconductor junctions. Of 
particular interest is the case where a semiconducting 
nano-wire with strong spin-orbit coupling is fabricated 
on top of a superconductor. The nano-wire can be tuned 
into a topological phase, supporting MBSs at its two 
ends, by applying a Zeeman field onto it and tuning 
the chemical potential into the Zeeman gap 0, H| . For 
this system the zero-bias conductance peak, recently ob- 
served in the tunneling current from a normal lead, would 
be considered as evidence for the presence of MBSs 0- 
[5J. Signatures of MFs in topological Josephson junctions 
@;0]j which consist of two superconducting leads coupled 
through a three-dimensional topological insulator, have 
also been reported . 

In a typical experimental setup for realizing MFs the 
host topological superconductor (TS) is placed in prox- 
imity to a top gate. An intriguing and experimentally rel- 
evant question is whether MFs, being inherently charge- 
neutral objects, can participate in the charge response of 
the system to the nearby gate. In this Letter we study 
the effect of a small time-dependent oscillatory compo- 
nent of the gate voltage on a mesoscopic system support- 
ing MBSs which is coupled to a single-electron bath. We 
show that, remarkably, the MBSs play a dominant role in 
the response to the oscillating gate voltage; and that, in 
particular, they give rise to a sub-gap dissipative charge 
response. More specifically, the low frequency charge re- 
sponse of the system is encapsulated in two quantities: 
an effective capacitance and a dissipative charge relax- 
ation resistance. We demonstrate that both are sensitive 
to the presence of MBSs, and contain unique signatures 



of the latter (see Fig. 2, Eq. (JT2J), Eq. (JT3J) and the discus- 
sion following them) . Retaining the Coulomb interaction 
is crucial for correctly identifying these quantities. We 
expect that these concepts will play an important role in 
future experiments attempting to measure signatures of 
MFs embedded in mesoscopic settings. 

The linear charge response to gate voltage oscillations 
for the setup of a quantum dot connected to a single lead 
has attracted considerable interest In the coher- 

ent tunneling regime the charge response defines a meso- 
scopic capacitance Co and a charge relaxation resistance 
R q according to 



C (l + iuC R q ), 



(1) 



where Q is the charge on the dot and V g is the gate volt- 
age. The capacitance C describes the static charging of 
the dot. It is generally different 1|| ljj from the geomet- 
rical capacitance C and depends strongly on the lead-dot 
tunneling width T. Energy dissipation is associated with 
a time scale CoRq, where Rq is known as the charge 
relaxation resistance. It assumes the quantized value 
Rq = h/(2e 2 ) which, unlike the classical limit, for a sin- 
gle channel junction does not depend on the transmission. 
This prediction remains valid also when effects of interac- 
tion are taken into consideration 11 , 12j ■ The quantum 



impurity models that have been studied so far to analyze 
the charge response include the Anderson model II], 13 1 
and the model for large metallic dots, described by the 
Ambegaokar-Eckern-Schon effective action [15 1. 

We now proceed to analyze the linear charge response 
in a different quantum impurity model, presented in 
Fig. 1, where a mesoscopic TS wire supporting two MBSs 
at its ends takes the role of the quantum dot. The 
wire is coupled to a top gate and, on one of its ends, to 
a normal metal (NM) lead. The potential of the normal 
metal lead is taken equal to zero (V = 0), while the po- 
tential of the gate V g (t) includes a small time-dependent 
perturbation. The Coulomb interaction on the dot is de- 
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FIG. 1. A schematic of a quantum RC circuit with Majorana 
modes (represented as two red dots) . A topological superconductor 
which hosts two MBSs at its ends is coupled via a tunnel contact 
(with hybridization tunneling width T) to a normal metal lead. The 
dissipative current through the tunneling contact is induced by the 
ac gate voltage V g (t). The charging energy strongly influences this 
current and the effective capacitance Cq. 



scribed via the non-zero charging energy E c — e 2 /2C 
where C is the geometrical capacitance. The presence 
of a charging energy introduces dynamics for the super- 
conductor phase which, combined with the coupling to 
the lead, induces phase fluctuations on the dot; there- 
fore, it influences the dissipation. More specifically, the 
two Majorana fcrmions at the ends of the wire 71 and 
72 (index 1 denotes the left side which is taken to inter- 
act with the metallic lead), with anti-commutation re- 
lations {7i,7j} = Sij, introduce the auxiliary fermion 
d = (71 + i72)/\/2- The number operator fid = d)d 
defines the occupation of this fermionic level. The two 
degenerate ground states with even or odd number of 
electrons are associated with different electron parities. 
This fixes the occupation of the fermionic level accord- 
ing to 2nd — 1 = (— l) n where n is the electron num- 
ber. Therefore, the MF response gets intermixed with the 
charge (or equivalently its dual phase) dynamics [l8j |. 

To calculate the charge response, we write the charge 
on the TS wire via the following formula 



Q(u) = e 2 K(uj)V a (io), 



(2) 



which is expressed in terms of the retarded response func- 
tion K(t - t') = i9(t - t') ([hd(t),n d (t')]), describing the 
charge fluctuations of the fermionic level at equilibrium 
(non-equilibrium response is also of interest [19| but is 
beyond the scope of this work). Here 9(t) is the step 
function: 6(t) = 1 if t > and is equal to zero if t < 0. 
Below we proceed similarly to the approach presented in 
[20j and construct a general Keldysh functional for the 
interacting problem using the dual phase representation 
for the Coulomb charging term. 

The Hamiltonian and the action. — We now derive the 
Keldysh generating functional in the presence of an ex- 
ternal source term, which later will be used to calculate 
correlation functions. Our starting point is the Hamil- 
tonian: The Majorana Coulomb box, lead and tunneling 



Hamiltonians are respectively [18|, l20( 

H c = E c (2N + n d - n ) , H L = J2k e kc\c kl 



H t = Act(0) (d + dt e - 2 ^) + h.c, 



(3) 



where no is tunable by the gate voltage V g and N is 
the number operator for Cooper-pairs, conjugate to the 
superconducting condensate phase 20. The lead Hamil- 
tonian H l is expressed in terms of quasiparticle ener- 
gies ek and electron operators Cfc,c£. The operators 
c(x),c'(x) are the Fourier transforms of the latter. By 
applying the Hubbard-Stratonovich transformation, the 
quadratic term in N is replaced by a gaussian form in 
</> with additional functional integration on the phases 
= (<ftii <fe) (where the phase 4> is half of the condensate 
phase) |2l|,|22j. The total Lagrangian and, specifically, 
the Majorana sector of the Lagrangian, both depend on 
these phases. This directly shows that the charge trans- 
port which is described by the phase variables strongly 
influences the Majorana dynamics [3, H(| ■ The action 
of the dot takes the form 



jdt 
Lm 



- + n <f> + Lm ) o- z , 
d{idt - X) d. 



(4) 
(5) 



Here x ~ <f> — a where a is a source field. In the above 
formulas and throughout this Letter the Pauli matrices 
o~o , <r x 7 o~ y , o~ z act within the Keldysh space. Accordingly 
all fields and sources have two components (denoted 1,2). 

It is convenient to gauge out the total field x (including 
the source components) from the action Lm by applying 
the transformation d — > de~ lx . In this gauge the tunnel- 
ing part of the action includes a phase factor 



Si 



A / dtc(0)cr z jie- zx + h.c, 



(6) 



while the pure Majorana part becomes independent of 



X : L M 



1 T 

2T 



[<7 Z G -1 ]7 where G is the Majorana Green 
function (GF). Following a Keldysh rotation this GF ac- 
quires the form (23| 



1 G A 



Each element of this matrix is itself a matrix in Majorana 
basis. Below we use only the 11-element of the latter 
matrices, being, 

Gfi(e) = 



Gfx(e) 



(e + i5) 2 - u 2 Q ' 

-wrtanh [6{e - vq) + S(e + vo)]. (7) 

Here in addition to the terms in Eq. ([5]) we included a 
coupling energy vq between the two MBSs. 

Integrating out the lead operators and performing the 
Keldysh rotations of basis (j Tl , -f')=(-f T o- z L~ 1 , Lcr z j), 



3 



where L = -^(fo — i (T y)i we obtain the effective action 
as function of the phase and source fields 



S eS ((f>,a) = i/ ' dtdt'^J'\G^ -2TS itj S ltj a x g x ]j' j , 



9x 



%t')e 



(8) 



( q R q K \ ■ 

Here 5 = ( g A I is the integrated over momentum 

matrix GF of the lead, while T = 2tt\ 2 N(0) stands for 
the tunneling width. All the fields arc presented in terms 
of their classical and quantum components as Xc = (xi + 

X2)/2 ;Xg = Xi- X2- 

The next step is the integration over the Majorana 
fcrmions. The resulting generating functional takes the 
form 



t>,a) 



cxp 



1. 



-Trln(G - 2T6 ilj Si ! ja x g x ] 



Integrating over the fields <f> we get the functional to a 
form which depends only on the source fields 



Z f (a) = J D((j) c ,(j) q )Z((j),a)( 



iS b 



Sb 



dt 



2E r 



(9) 



(10) 



We are looking for a perturbative expression for InZ/ to 
4th order in tunneling: 



HZ f {a)]=Hzf\a)]+HZf{a)l (11) 

where the two terms in the right hand side of Eq. (fTTj) 
correspond to second and fourth order contributions re- 
spectively (see the supplement for details). 

Effective capacitance. — The effective capacitance Co 
does not coincide with the geometrical capacitance C and 
is expressed to second order in tunneling through the 
response function Cq = e 2 RcA'(oj = 0) in the the 
static limit. Generally the real part of K(t,t') takes the 
form 



(2), 



RcK(tt') 



2 \Sa c (t)6a q (t') Sa q (t)5a c {t') 



\n\zf\a)]. 



Lengthy but straightforward calculations (see Supple- 
ment) lead to the following result in the limit of zero 
temperature (T — > 0) where we concentrate on the near- 
est two peaks 



C = e 2 RcA(u; = 0) 



CT x - 1 



,(12) 



where n v = 1 + 2tiqv. 

Thus Co generally depends on the interaction energy 
vq between the two MBSs. For ^ 0, the degeneracy 
points are no longer equally spaced (with a spacing equal 
to e) but they bunch in groups of twos due to alternate 
spacings e(l±^o/A c ) (see Fig. 2 where the entire series of 



peaks is presented and not only the two resonant ones). 
See also [HHHI for a quantum Hall analog of this effect. 
Charge spectroscopy, which is realized by measuring the 
effective capacitance, can therefore display revealing sig- 
natures of the MBSs. 




■ H 



FIG. 2. The effective capacitance Co is displayed schematically as 
function of the background charge no for the Majorana Coulomb 
box problem with spacings that approach e when the Majorana 
interaction — > (solid blue line) versus a regular superconductor 
with equal spacings 2e (dashed black line), assuming A 3> E c . 

Charge relaxation resistance. — The imaginary part of 
density-density correlation function Imif(tt') is obtained 
by the variation of In Zf with respect to the source fields. 
In the limit of small frequency (uj — > 0) the relevant 
contribution to ImK comes from the fourth order term 

,(2) 



HZf\a)] 



ln[Zy>] of Eq. (HU) 

I f s 2 s 2 

m X{U ' - 2 \Sa c (t)Sa q (t') 5a q (t)Sa c (t>) 

To calculate the imaginary part of retarded correlation 
function we proceed similarly to the case of the real part. 
Using the second term of Eq. (fTl"]) the frequency depen- 
dent imaginary part of response function ImAT(w) (see 
supplement) acquires a form 



ImA'(w) = £ £ 

w' — ± 



de 1 de 3 0(i'ei)0(v'e 3 )G-(e 1 )G-(e 3 ) 



where G_ = G{\ — . Completing the integration over 
energy variables we obtain the main result of this work 



ImA» = ^ 



8tt£V 



,i n i n i> 



(13) 



where n v — n v + i>o/E c and we restore the units h = 2n. 
Using the last formula we see that Im A(w) satisfies the 
Korringa-Shiba relation [26| 



ImA'(w) 



^[AeA(O)] 2 . 



Thus we demonstrated the presence of a dissipativc 
charge relaxation resistance in response to the oscillating 
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gate voltage, associated with a unitary resistance value the simple quadratic action 
R q = h/2e 2 . This result deviates from its expected value 



when the interaction on the dot is taken to zero (E c = 0). 
The value for the charge relaxation resistance of the Ma- 
jorana Coulomb box is equal to that of a local impu- 
rity level which is described by the Anderson Hamilto- 
nian Superficially the Majorana Coulomb box in the 
presence of strong Coulomb blockade mimics the physics 
of a small quantum dot coupled to a metallic reservoir. 
The Majorana box discussed here is nevertheless different 
in several important aspects: (i) the superconducting gap 
is the largest energy parameter of the system (A 3> TV) . 
Therefore the critical Fermi liquid hamiltonian (see |14| ) 
is not reached and the Kondo effect is excluded; (ii) in 
addition to regular tunneling, there is also an anomalous 
(Andreev) term in the tunneling hamiltonian ([3]) which 
is absent in the standard tunneling hamiltonian of the 
Anderson model. The first one describes the transfer of 
an electron from the d-level to the normal lead, and the 
second one the annihilation of a Cooper pair inside the 
island accompanied by the creation of two electrons, one 
in the d-level and one in a normal electrode; and, (iii) 
the occupation fid of the d-level depends on the fermion 
parity which imposes a constraint on the Hilbert space 
([HI]). We dismantled this constraint (see [lil) by link- 
ing the occupation of the Majorana bound state to the 
parity of the superconducting condensate. 

The principal conclusion that follows from the differ- 
ences between the Majorana Coulomb box described here 
and the Anderson model is that for the Majorana case 
the strong interactions influence the conductance (damp- 
ing the Andreev tunneling) and, unlike the case without 
interactions, result in the same value for charge relax- 
ation resistance as that of a quantum dot with a single 
energy level. This is in stark contrast to the case of the 
Anderson model where the charge relaxation value is the 
same independently of interactions. 

Resonant tunneling model. — Equations (|12p and (|13|) 
have been derived using the Keldysh formalism. It is 
useful to compare them with an effective resonant tun- 
neling model obtained in [l8|| describing the Majorana 
box in the strong Coulomb blockade regime. The effec- 
tive Hamiltonian focuses on two charge states that differ 
by one electron [ljj ]. For our setup (Fig. 1) it takes the 
form 



H = Hi 



+ eo (V/-^) +Act(0)/ + h.c. 



(14) 



where /' (/) are the creation (annihilation) operator for 
a Dirac fermion occupying the resonant level, eo is the 
energy difference between the levels of resonant model 
and depends on the gate voltage. The density-density 
correlation function K(t) is expressed in terms of the 
occupation n/ = f^f of the resonant level. After we 
integrate out the reservoir the problem is described by 



dtfGj\f. 



(15) 



Here the matrix form of the GF G/ is the same as the one 
for g in Eq. ([8]) and its entries are: Gj (e) = Gy * = (e — 
eo + iT/2)- 1 , Gf(e) = i rtanh(e/2T)/[(e-e ) 2 + r 2 /4]- 1 . 
With the help of these Green functions we can immedi- 
ately find the dynamical charge response and effective 
capacitance (see details in the supplement). In the low 
frequency limit we obtain 



ReK(uj -> 0) 



r 



i 



2ir e 2 
hcu , 



r 2 /4' 



lmK((j) = ^[Re/-iT(0)] 2 . 



(16) 
(17) 



By comparing the results of our direct expansion in tun- 
neling, Eqs. (fr2|) - ([T5|) . with the results obtained using 
the resonant model, Eqs. (|TrI|) - p7|) . we can find the gate 
voltage dependence of eo- To second order in tunneling 
we get 



eo 



E, 



E 



2(n v + u Q /E c f 



-1/2 



(18) 



Hence the results derived using the effective resonant 
model described in 18| are fully consistent with our di- 
rect perturbation theory approach once the effective pa- 
rameters are correctly identified. 

Conclusions. — In this work we used the Keldysh tech- 
nique to investigate the dynamical charge response in a 
system formed by a short topological superconducting 
mesoscopic wire with two Majorana fermions hosted at 
the ends of the wire. The wire is coupled on one side to a 
normal lead and low frequency gate voltage is applied to 
the system. We calculated the effective capacitance and 
charge relaxation resistance and derived a dissipative re- 
laxation resistance associated with a unitary conductance 
value R q = h/(2e 2 ). This is in agreement with Korringa- 
Shiba formula where the charge relaxation resistance is 
equal to R q . We demonstrated that this value is also con- 
sistent with the effective model for the strong Coulomb 
blockade limit described in [l|| and identified its effec- 
tive parameters. The capacitance Go generally depends 
on the interaction between the two Majorana states (see 
Fig. 2). Both these effects (dissipation and capacitance) 
can help identify Majorana fermions. 
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SUPPLEMENTARY INFORMATION 
Perturbation Theory 

The second and fourth order terms in Eq. (11) of the 
main text consist also of a functional integration over the 
phase (see Eq. (9)) 



\n\zf\a)] = -r(TrO(tt)), 



(SI) 



where G u (i/,ei) - [G& + u(G& - G£)] fl , g K (e) = 
— i tanh(e/2T). Finally, following integration over en- 
ergy variables and taking the limit T — > we arrive at 
the main result for the effective capacitance presented in 
Eqs. (12)-(13) of the main text. 

Charge relaxation resistance: Charge relaxation resis- 
tance can be extracted from the fourth order term of 
Eq. (11). We start our calculations by writing the ex- 
plicit form of the second term in (|S2I) : 



i((TrO(tt)) 2 - ((TrO(tt)) 2 )) + Tr [flfaMht)] = ±f ^UMtdt^dts 



(Tr 



(S2) 
(S3) 



Here the trace includes an integration over the time vari- 
ables. 

Effective capacitance: To find the effective capacitance 
we need to calculate the real part of the correlation func- 
tion K(u> — > 0). Explicitly the second order contribu- 
tion Eq. (|S2j) (prior to the functional integration over 
the phase) acquires a form 

Trfi(tt) = i / d i ^ 1 ^e*[^(ti)-xc(t)+i(x 9 (ti)— x«(*)-')/2] 

v.v' 

xTr[Gn(i£iK(l + va x )g{t x t){\ + i/<r x )].(S4) 

Here u, v 1 = ±1. We notice the relation Xe,g(ii) ± 

Xc,,(*) = 0c,,(ti)±0c,,(t)- Jo 1 dt'a Ctq (t')± f*dt'a Ctq (t'). 
Variation with respect to classical source field a c and in- 
tegration over <f> c define a functional 5-function which sets 
the value of the quantum field cj> q : 



(S5) 



<j> q {t') = 2E c [6(t 1 - t')(h - t') - 0(t - t')(t - t')}. (S6) 

Variation with respect to the quantum source defines 
KeK (using Eq. (|S6[1 ). The principal contribution to the 
variation of In Zf comes from v 1 = — v and is equal to 

5 2 Hzf ) ] a ^ Q {t } ) = ^/Z j dtdtx(6(t) - 0(-ti)) 

X e -ivE u (t-ti) p 

x[e(ti-t f ) + e{t-t f )]e(tfUsr) 

where E v — E c (l — 2n$v) and we denote the trace 
term in Eq. (|S4|) as F. We Fourier transform according 
to J dtfexp[iujtf]..., and in the limit of zero frequency 
the term in the square brackets is integrated to a form 
tid(ti) + t0(t). Direct integrations over t, t% now result 
in: 

ReK = _E f 'ti^i v ei)[gK ~ v{gR ~ gA)]t2 



{E v + v{ex - e 2 )) 3 



where 



ni = cxp {-i[ X c(t) + Xc(h) - Xc(*a) - Xc(t 3 )}} : 



and 

n 2 = 



cxp 



Here the integration is over all time variables and sum- 
mation refers to all numbers = ±1 (i = 1, 2, 3, 4). The 
phases x are the sum of charge and source fields, and 

*(„ t ) = Tr[G , n(ei)cr a; (l + v 1 <j x )g(e 2 ){l + v 2 a x ) 

Gn(e 3 )a x (l + v 3 <r x )g(e 4 )(l + i/ 4 cr x )]. (S9) 

The above ^> is analogous to the function F introduced 
before. Performing variations with respect to the quan- 
tum and classical source fields we arrive at the following 
expression for ImK (ui) 

ImK(u) = Llm J J2[9(t) + 9(h) - 9(t 2 ) - 9(t 3 )] 

x* (i/ . ) e^[ £,+2n °( t+tl -* 2 - t3 ^r(w), (S10) 

where the terms in square brackets originate from the 
variation with respect to the classical source, while the 
function 



r(w) = -±-(-v x e 1 ^ 
2iuj 



v 2 e 



icjti 



v 3 e 



appears as a result of the variation with respect to the 
quantum source a q and is dependent on the path fol- 
lowed by the electron. The functional integration on the 
classical component of the charge field <p c with the action 
given in Eq. (10), defines the function D in the exponent 
of Eq. (|S10|) . To calculate D we use the solution for the 
quantum field <pq which is found from an equation similar 
to Eq. (|S5j) and results from integration over <f> c , 



= - 



+ 5{t 2 - t') - S(h - t') + 



(S8) 



2E C 

6(t 3 -t')-6(t-t'), 
<f> q (t') = 2E c [6{t 2 - t')(t 2 - t') - 9(ti - t')(h 
9{h ~ t')(t 3 - t') - 6{t - t')(t - *')]• 



t') 



(Sll) 
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Finally, replacing the charge quantum phase by its so- 
lution (fSTTj) we find the function D appearing in the 
exponent of (|S10[) . For a general path followed by the 
electron, D takes a form 

D=(t 2 - ti)M(*i - h) - v 2 6{t 2 - t x ) - e 2 ] + 
(*s - h)WiO(t 3 - t 2 ) - v 3 0{t 2 - t 3 )] + 
(t 2 - t)[v!e(t - t 2 ) - v 4 8(t 2 -t) + ei] + 
(ti - h)[v 2 6(h - tx) - v 3 6{h - t 3 ) - e 3 ] + 
(t - h)[v 2 6{t - h) - ViQiti - t)] + 
(t 3 - t)[v 3 6{t - t 3 ) - v A 6{t 3 - 1) - e 4 ]. 

In the definition of D we have included the energy vari- 
ables (in units of E c ) which come from the Fourier images 
of the Green functions [see in Eq. ([S9|l ]. Thus D 

defines the path in time which the electron follows. We 
identify each path by four numbers as ( v\v 2 v 3 V/±). The 
electron-hole excitations in the normal metal lead pro- 
vide the main contribution to the dissipation. This fixes 
4 relevant path configurations: (a)(— 111 — 1), (6)(1 — 1 — 
11), (c)(-l - 1 - 1 - 1), (cf)(mi). For these paths the 
trace in Eq. (|S9[) (taking the relevant part which con- 
tributes to ImK) acquires the following form at T — > 



40(± e2 )0(TC4)[Gfi(ei)Gfi(e 3 ) - G + (ei)G + ( C3 )] 



(a),(h) 

G + (ex)G+(e 3 )]. 



(S12) 



Here and below G± = G* ± Gft. In Eq. (fST2|) we 
kept only terms which have poles in the particle-hole po- 
larization (e 2 — £4 ± oj)^ 1 . For the four relevant paths 
(a), (b), (c), (<i) we also calculate the function D and find 
the corresponding contribution to ImK by performing 
the time integrations in Eq. (|S10j) . The time integrations 
fix the signs of energy variables and result in a factor 



E 

vv'=±l 



(S13) 



Combining the GF with Eq. (|S13|) and collecting all con- 
tributions we arrive at Eq. (14) for ImA'(oj) of the main 
text. 



Effective resonant model 

An effective resonant model describing the physics of 
the Majorana box in the strong Coulomb blockade regime 
is given by the Hamiltonian in Eq. (14). After integra- 
tions on reservoir variables the problem is reduced to 
the pure quadratic action Eq. (15). This action is taken 
on the Keldysh space (either directly or rotated by the 
matrix L of the main text). In direct space the correla- 
tion function is presented by time ordering products of 
n/ = /'/ operators 



K R (t - = i{T[n ( p{t)n\ L} (t') - nf>(t)nf> (*')), 



(i) 



,( 2 )/ 



ReK R (t - t>) = \{T[nf{t)nf(t') - nf \t)nf '(*')> , 



,(2),.s (2), 



lmK R (t-t') = ^{T[nf{t)n^{t') 



where the numbers in brackets denote the Keldysh in- 
dices on the time contour. Using Wick theorem we can 
contract the /-operators and express the above values in 
terms of the associated Green functions G/ 



ReK R (uj -> 0) 



de ft 
— tanh — - 
2tt \2T 



V ^ 4 7 2tt 



tanh 



2T 



G/_(e)G /+ (e), 



— tanh ( — ) 
27V 



where the form of the Green functions G/± = G R ± Gj 
and Gf is presented in the main text. Direct integrations 
over the energy variables result in Eqs. (17), (18) of the 
main text. 



